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Magnetic field dissipation in neutron star crusts: 
from magnetars to isolated neutron stars. 
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ABSTRACT 

Context. We study the non-linear evolution of magnetic fields in neutron star crusts with special attention to the influence of the Hall 
drift. 

Aims. Our goal is to understand the conditions for fast dissipation due to the Hall term in the induction equation. We study the 
interplay of Ohmic dissipation and Hall drift in order to find a timescale for the overall crustal field decay. 

Methods. We solve numerically the Hall induction equation by means of a hybrid method (spectral in angles but finite differences 
in the radial coordinate). The microphysical input consists of the most modern available crustal equation of state, composition and 
electrical conductivities. 

Results. We present the first long term simulations of the non-linear magnetic field evolution in realistic neutron star crusts with 
a stratified electron number density and temperature dependent conductivity. We show that Hall drift influenced Ohmic dissipation 
takes place in neutron star crusts on a timescale of 10 6 years. When the initial magnetic field has magnetar strength, the fast Hall drift 
results in an initial rapid dissipation stage that lasts ~ 10 4 years. The interplay of the Hall drift with the temporal variation and spatial 
gradient of conductivity tends to favor the displacement of toroidal fields toward the inner crust, where stable configurations can last 
for ~ 10 6 years. We show that the thermally emitting isolated neutron stars, as the Magnificent Seven, are very likely descendants of 
neutron stars born as magnetars. 
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1. Introduction 

Since the early days of neutron star (NS) research, when a 
first estimate of the c haracteristic deca y time of its magnetic 
field was performed (iBaym et al. 1969I) . our models of mag- 
netic field evolution in isolated NSs has become more and 
more complex. However, a complete theoretical model that ex- 
plains all observational facts satisfactorily does not yet exist. 
There is quite general agreement that the magnetic field in ultra- 
magne tized NSs (magnetars) decays on timescales of 1 3 - 1 5 
years (IWoods & Thompson 20061; lHarding & Lai 20061) . Very 
recently, it has been argued (Po ns et al. 20071) that field decay 
may be an effective heating mechanism also in isolated NSs 
which are somewhat less magnetized than standard magnetars 
(AXPs and SGRs). The observed correlation between surface 
temperatures and magnetic field strength is an evidence for effi- 
cient dissipation, and hence Joule heating, near the surface lay- 
ers of NSs. In addition, the evidence for braking indexes signif- 
icantly larger than 3, which are inferred for an ensemble of a 
dozen of pulsars in an active age of a few 10 5 years, indicate that 
epochs of rapid decay in relatively young NSs occur (Johnston 
& Galloway 1999, Geppert & Rheinhardt 2002). 

Another evidence for rapid magnetic field evolution in the 
subsurface crustal layers would be the observation of small 
structures (significantly smaller than the dipolar mode) close to 
the surface. These small scale structures could not be survivors 
from the magnetic field structure acquired by the NS at birth; 
their small scales, in combination with the low electric conduc- 
tivity in the outer crust, would have erased such field structures 
on timescales shorter than 1 Myr. The existence of small scale 
structures is, however, a necessary condition for the pulsar mech- 



anism to work (Ruderman & Sutherland 1975). In addition, the 
growing understanding of the physics of the drifting sub-pulses 
supports the idea that small scale crustal field structures must be 
present (Gil et al. 2003, 2006). The existence of strong crustal 
fields in NSs as old as about 1 Myr is also necessary to cre- 
ate the non-uniform surface temperature distribution observed 
e.g. in the "Magnificent Seven" (see Pons et al. 2002, Schwope 
et al. 2005; Haberl et al. 2006 for the observational evidence 
and Geppert et al. 2004, 2006; Perez- Azorin et al. 2006a, 2006b 
for the theoretical interpretation). All this phenomena can only 
be understood if there is a crustal magnetic field with a strong 
toroidal component, and with more structure than the simple 
dipolar model. 

On the other hand, population synthesis studies suggest 
that old pulsars show no significant magnetic field decay over 
their life time, i.e. the decay time must be larger than 10 Myr 
(Hartman et al. 1997; Regimbau & de Freitas Pacheco 2001), al- 
though the opposite conclusion has also been claimed (Gonthier 
et al. 2004). These, at first glance, contradictory facts can be sat- 
isfactorily resolved by the (quite natural) assumption that the NS 
magnetic field is maintained by two current systems. Long liv- 
ing currents in the superconducting core support the large scale 
dipolar field and are responsible for the spin down of old pulsars. 
Currents in the crust support the short living part of the field. It 
decays on a timescale of 10 5 — 10 7 years, depending on the con- 
ductivity, thickness of the crust, and strength and structure of the 
initial field. Estimates of how fast a core anchored field could be 
expelled and subsequently dissipated in the crust result in char- 
acteristic timescales exceeding 100 Myr (Konenkov & Geppert 
2001). The dipolar component of the crust is superimposed, out- 
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side the NS, to the core component. Depending on their relative 
strengths, a rapid decay of the crustal dipolar field may or may 
not have observable influence on the pulsar spin down behaviour. 

Besides the Ohmic diffusion, which will proceed fast in 
the outermost low-density crustal regions and during the early 
hot phase of a NS's life, the only process that can change 
the crustal field structure, both quantitatively and qualitatively, 
is the Hall drift. Many authors studied during the last two 
decades the effects of Hall drift onto the evolution of mag- 
netic fields in isolated NSs (Haensel et al. 1990; Goldreich 
& Reisenegger 1992; Muslimov 1994; Naito & Kojima 1994; 
Urpin & Shalybkov 1995; Shalybkov & Urpin 1997; Vainshtein 
et al. 2000; Rheinhardt & Geppert 2002; Geppert & Rheinhardt 
2002; Hollerbach & Ruediger 2002; Geppert et al. 2003; 
Rheinhardt et al. 2004; Cumming et al. 2004). By use of the anal- 
ogy of the Hall induction equation with the vorticity equation 
of an incompressible liquid, Goldreich & Reisenegger (1992) 
developed the idea of the Hall cascade. It transfers magnetic 
energy from the largest to smaller scales until a critical scale 
length below which Ohmic decay becomes dominant is reached. 
Vainshtein et al. (2000) considered for the first time the Hall drift 
in a stratified NS crust. In the presence of a density profile, Hall 
currents are able to create current sheets, which are places where 
very efficient dissipation occurs. This effect is dramatic if current 
sheets are located just below the surface, where the conductiv- 
ity is lower, but even if the drift is directed towards the highly 
conductive inner crust, the small scale of the locally intense field 
causes a significantly faster dissipation of magnetic energy than 
the purely Ohmic diffusion estimate. In summary, this MHD- 
like process (considering the crust as a one-component plasma) 
is by itself energy conserving. However, it affects the magnetic 
field by two inherent tendencies: the creation of small scale field 
structures by transferring magnetic energy from the initial large 
scale field and, the drift of growing structures towards a region 
where current sheets are created. 

Another interesting effect of the Hall drift was proposed by 
Rheinhardt & Geppert (2002). They show, by means of a sta- 
bility analysis of the linearized Hall induction equation, that the 
transfer of energy from the large scale (background) field to the 
smaller scale modes may proceed in a non-local way in the mo- 
mentum space, resulting in a Hall instability. Cumming et al. 
(2004) showed that the growth rate of this instability depends on 
the shear in the velocity of the electrons whose current supports 
the background field. As they pointed out it is unclear whether 
the Hall instability is relevant for the field evolution in the crust, 
since the Hall cascade may proceed sufficiently fast to fill in the 
intermediate scales between the large scale initial field and the 
unstable growing small scale ones. 

In the present study we intend to consider the field evolu- 
tion as realistic as possible. Thus, we will solve numerically 
the non-linear Hall induction equation in the crustal region, us- 
ing a crustal density and conductivity profile according to the 
state-of-the-art microphysical input. We start with initial field 
configurations which consist both of poloidal and toroidal field 
components that fulfill boundary conditions at the superconduct- 
ing core and at the surface. We will follow the cooling history 
of a NS, starting with an initially hot NS where the magneti- 
zation parameter is relatively small and the field evolution in 
the crust is almost completely determined by Ohmic diffusion. 
As the NS cools down, the magnetization parameter increases, 
thereby gradually enhancing the relative importance of the Hall 
drift. Our aim is to reach a better understanding of the effects that 
Hall drift may have on the crustal field evolution and its conse- 
quences for the rotational and thermal evolution of isolated NSs. 



The paper is organized as follows. In Sect. 2 we present the 
Hall induction equation where the magnetic field is represented 
in terms of its poloidal and toroidal components which are ex- 
panded in a series of spherical harmonics. Sect. 3 is devoted to 
purely toroidal fields. For this case we show that the induction 
equation can be written in a form similar to the Burgers equation. 
This allows us to discuss the importance of the field-dependent 
drift in a clear way. In Sect. 4 we detail the input microphysics, 
the NS model, and the initial conditions. In Sect. 5 we present the 
results for different initial field configurations. Finally, in Sect. 
6, we discuss our main conclusions and give an outlook to future 
improvements. 

2. Basic Equations and Formalism 

In the crystallized crust of NSs, where convective motions of the 
conductive material play no role, the evolution of the magnetic 
field is governed by the Hall-induction equation 

dB c 1 1 1 \ 

— = -— Vxl-{VxB + ^r[(VxB)x6]}l, (1) 

where b is the unit vector in the direction of the magnetic field 
b = B/B, with B being the magnetic field strength, t is the relax- 
ation time of the electrons and u> B = eB/m* e c is the electron cy- 
clotron frequency. Here, <x = e 2 n e T/m* is the electrical conduc- 
tivity parallel to the magnetic field , and n e the electron number 
density. The Hall drift term (the term proportional to wjt) at the 
r.h.s. of Eq. (Q]i is a consequence of the Lorentz force acting on 
the electrons. The tensor components of the electric conductiv- 
ity are derived in the relaxation time approximation (Yakovlev & 
Shalybkov 1991). If the magnetization parameter wgT exceeds 
unity significantly, the Hall drift dominates, which results in a 
very different field evolution from the purely Ohmic case. A 
large ma gnetization parameter, typically » 1000 during som e 
stages dGeppert & Rheinhardt 20021: iPerez-Azorin et al. 2 006). 
strongly suppresses the electric conductivity perpendicular to the 
magnetic field. 

The up to now limited success to solve Eq. (Q]i has 
been restricted to magnetiza tion parameter not exceeding 200 
dHollerbach & Rudiger 2002) Q. These numerical limitations 
guided us to use an alternative approach, since fully spectral 
codes have systematically unsurmountable problems to deal with 
structures where discontinuities or very large gradients of the 
variables appear. This will happen, for example, if current sheets 
develop as a consequence of the Hall drift. The next section will 
describe in detail the reason why spectral methods may not be 
well suited to deal with this problem. We have decided to em- 
ploy a semi-spectral method which describes the angular part 
of the field by spherical harmonics but uses a spatial difference 
scheme in radial direction. 

The formalism that uses the representation of fields by their 
poloidal and toroidal parts and their expansion in a series of 
spherical harmonics was developed by Radler (1973). Radler's 
formalism allows to transform the (vector) induction equation 
into two scalar equations. In this formalism, the magnetic field 
is decomposed in poloidal and toroidal components: 

B = B po ] + B toI . (2) 

Hereafter, we follow the notation of Geppert & Wiebicke (1991). 
The two components are described by two functions, (£>(/-, 8, (p, f) 

1 Their parameter R B is equivalent to our lo b t calculated for a fixed, 
initial magnetic field strength So. 
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and *P(r, 6, ip, t), where r, 6, and ip are the usual spherical coordi- 
nates. Explicitly, the components of the field are given by 



Bpoi = V x A , 
Bpoi = = -r V 2 <5 + V 
B tnl = -rxVF. 



A = -r x V<1> 

d(r®)) 



(3) 



This formulation has the advantage that it automatically fulfills 
the divergence condition (V ■ B — 0). Inserting the expressions of 
Eq. ([3]) into Eq. (Q]) we arrive at two partial differential equations 
describing the evolution of the poloidal and of the toroidal part 
of the magnetic field. 



<90> 


1 


~di ~~ 


& 




1 


~di~ ~~ 


& 



1 dlogo-dfW) 
r dr dr 



+ C . 



(4) 



Here, & = 4ncr/c 2 and D and C stand for the nonlinear terms de- 
scribing the Hall drift and coupling both components (poloidal 
and toroidal) of the magnetic field. We have assumed that the 
conductivity depends only on the radial coordinate, and is inde- 
pendent of the magnetic field strength. 

After expanding the functions O and *P in a series of spheri- 
cal harmonics 



$ = -Yj^nm(r,t)Y nm {e,4>) , 

n,m 

x ¥=~Y J x i'nm(r,t)Y lm (6,cf>), 



(5) 



where n — 1 , . . . , « max and m = -n, . . . , +n. The vector potential, 
as well as the poloidal and toroidal parts of the magnetic field, 
can be written as 



A - - 



\ n,m 



dY 



~8~r ~d6 



eg , 



1 \ 1 dY nm 



(6) 



while the corresponding components of the current density, 
given by J = x B, is 



An 1 1 

—J = — 2j + V^nmYnm e, +- ^ 



72 



n(« + 1) 



dr 2 



11,111 

dY 



" A nm u± nn 

~dr dJ 



■ eg 



(7) 



Finally, by inserting the expansions of Eq. (|6]l into Eq. ([]]), we 
arrive to an infinite set of partial differential equations: 



dt 
dt 



J_ / d 2 Q„ m »(»+!) 
a" \ <9r 2 r 2 

d" \ dr 2 



+ D„ 



»(» + l) Uf 1 d&ffV nm 
r z cr or or 



(8) 



In this paper we will restrict ourselves to axially symmetric field 
configurations, i.e., the index m — and we will drop it hence- 
forth. Following the derivation of Geppert & Wiebicke (1991) 
the nonlinear coupling terms are 



(9) 



°~ r f~p \" r dr 



— y 

or 2 f-* 

k,k' 



r(3) 



_3 



where 



(10) 



(11) 



Above, f = cj b t/B, and / (2) J (3) are expressions that contain the 
Clebsch-Gordan coefficients which reflect the coupling proper- 
ties of the field modes with different multipolarity (see Eq. (60) 
in Geppert & Wiebicke (1991)). 

Using the orthonormality properties of the spherical harmon- 
ics, the volume integrated magnetic energy can be calculated as 



dV B 



1 J dr £>(n + l) 



.(12) 



As a criterion for the magnetic energy conservation we will 
check that the equality 



1 d 
8n dt 



J B 2 dV = - j J — dv > 



(13) 



is satisfied during the evolution within a certain tolerance (typi- 
cally less than 10~ 3 ). The integral on the l.h.s is evaluated over 
the volume occupied by the magnetic field (including the vac- 
uum exterior region), while the integral on the r.h.s. is calculated 
over the region where currents exist (the crust). 

2.1. Outer Boundary Conditions 
2.1 .1 . Outside Vacuum 

Since we consider in this study realistic crusts with finite elec- 
tric conductivity, surface currents are excluded. This means that 
we require all components of the magnetic field to be continu- 
ous across the NS surface, i.e. that the scalar fields <S>„ and 
and the derivative ^ are continuous through the outer boundary 
(see also Radler 1973). 

The external vacuum solution of a NS magnetic field is de- 
termined byVxB = 0, V-B = and the boundary conditions. 
From Eq. (0, the vanishing curl condition leads to 



<9 2 <D„ n(n + 1) 



dr 2 



<S>„ 



dY n 

-^e = at r>R, (14) 



which can be expressed as A<D = 0, where A is the Laplacian. 
The index m = in Y„ m has been omitted for simplicity. The 
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only physical solution of this equation is <!>„ = a„r ". Therefore, 
the requirement of continuity across the surface results in 



—1 = — $ at r = R. 

dr R 



(15) 



Because the poloidal current must also vanish in vacuum, we 
can derive another general boundary condition 



z 



3(sin0ff) 



rsinO 



do 



+ -■ 

r 



dr 



-eg 



= for r >R .(16) 



The existence of surface currents may affect the condition on the 
^-component of the poloidal current density, but there cannot 
exist radial currents penetrating into the vacuum. Thus the outer 
boundary condition for the toroidal field is simply 



= at r = R. 



(17) 



In general, the poloidal tangential ^-component of the cur- 
rent density has to vanish at the surface only if the electric con- 
ductivity vanishes there. Then, according to Ohm's law, any fi- 
nite tangential current density would cause an infinite tangen- 
tial electrical field which is in contradiction to the energy con- 
servation guaranteed by Maxwell's equations. Therefore, in the 
case of cr = at r = R, any solution of the induction equa- 
tion which fulfills the boundary condition W„ = at r — R 
will be characterized by vanishing tangential surface currents, 
i.e. ¥„ = = at r = R. 



2.2. Inner Boundary Conditions 

The inner boundary conditions are determined by the transition 
from normal to superconducting matter at the crust-core inter- 
face r = Rj. The Meissner-Ochsenfeld effect demands that the 
normal component of the magnetic field has to vanish at r = 
Continuity of the tangential component of the electric field to- 
gether with Ohm's law enforces that component to vanish at 
r = Rj because otherwise the infinite electric conductivity would 
cause infinite tangential current densities, finally destroying the 
superconducting state. 

For a spherically symmetric NS, the normal component of 
the magnetic field is its r-component as given by Eq. ©. Thus, 
the inner boundary condition for the poloidal field is 



<!>„ = at r = Rj 



(18) 



The tangential component of the electric field consists of 0— 
and ^-components. For the (/^-component we find, after some 
algebra: 



-Z 



»(» + !) „, dgy 8Y n , 
r 3 " dr " dO 

<9 2 <D„ n(n+l) 



n'(n' + l)ff¥„ ^ SY^ 
r 3 dr "' dO "' 



roj B T 



Z 



dr 2 



dO 



This electric field component has to vanish at the surface of 
the superconducting core. Since O,, = at r — R, , this condition 
reads 



s 



n(n + iyV n — — Y„ — 

dr dO 



-E 



d 2 <S>„ dY n 
dr 1 dO 



= 0.(20) 



1000.0 



100.0 



10.0 



1000.0 




Fig. 1. Radial profiles of electrical conductivity cr (solid lines), 
and the magnetization parameter normalized to a fiducial mag- 
netic field ofBn = 10° G (dashed lines), for three different 
temperatures. The crust is assumed to be isothermal. 

This condition is obviously not suited to find an inner bound- 
ary condition for the toroidal field. Therefore, we have to con- 
sider the ^-component of the electric field which consists of 
three contributions: 



(V x B tor ) x B t0l + (V x B pol ) x Bpo, + —(V x B tor ) 

(JJ B T 



21) 



Using again the expressions of Eq. © and taking again into 
account that <5„(r = Rj) = 0, the term (V x B po i) x B po i vanishes 
and it remains from the condition Eg = at the inner boundary 
that 



u>bt t—i dr dO r L 



dY„ 



(22) 



8Y] 



Multiplying both sides with -^4- and integrating over the solid 
angle the orthonormality of the spherical harmonics returns 



1 &¥„ 
cobT dr 



k,k' 



(23) 



Hollerbach & Rudiger (2002) applied the above boundary con- 
dition in the limit of lobt — » oo, avoiding thereby the difficulties 
the non-linearity will cause. We use the general form of Eq. j23l , 
since during the cooling process of the NS, there certainly will 
be periods during which it is not justified to neglect the dissipa- 
tive term. 



3. Comments on the evolution of purely toroidal 
fields 



From Eqs. ©, (0, and ( TTOb it is readily seen that an initial purely 
(19) toroidal field (<D = at t — 0) will not develop a poloidal part 
and does remain purely toroidal. For clarity, we will use in this 
section the (,0-component of the magnetic field (B^) as a variable. 
In axial symmetry, it is related to *T through 



_5 V F 



do eie 



(24) 



In order to make evident the effect of the Hall drift we will 
use cylindrical coordinates (R, ip, z). In the Appendix A we have 
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written the general form of the induction equation in cylindri- 
cal coordinates, with a decomposition of the magnetic field as 
used by Hollerbach & Riidiger (2002). When considering the 
case of purely toroidal fields, neglecting for simplicity the dissi- 
pative term (limit of strong magnetization), and assuming con- 
stant electron density n e , the induction equation is reduced to a 
single evolution equation for the variable B^. 



dB^ t 

— 1 = -- e ■ V x [(V x B tor ) x B tor ] 
at cr 



2B^ dBy, 



4nen e R dz 



(25) 



where we have used that t/cf — clAnen e . 

In this form, the Hall-dominated induction equation has the 
form of the inviscid Burgers equation with a wave velocity that 
depends on the variable B^, and on the coordinate R. This has 
been pointed out by Reisennegger et al. (2005) who also con- 
cluded that it leads to the formation of current sheets. In this 
way, the formulation illustrates very clearly that the originally 
(in the limit of weak field) purely parabolic diffusion equa- 
tion changes its character to hyperbolic when the Hall term 
dominates. The analogy between the induction equation and 
the Burgers equation was discussed before by Vainshtein et al. 
(2000) but there are a few more qualitative differences we should 
notice. Vainshtein et al. (2000) derived the analogy with the 
Burgers equation in Cartesian coordinates, assuming no depen- 
dence of the field on the z-coordinate. In this particular case, it 
is necessary to consider a stratified medium in order to have the 
Hall drift leading to the formation of current sheets. Our result 
is more general: even in the uniform density case, the geometry 
of a NS crust (axially symmetric toroidal field in a conducting 
spherical shell) leads to a Burgers-like equation that does not ad- 
mit stationary solutions. Depending on the sign of the gradient, 
the field will drift vertically in one or the other direction, lead- 
ing to the formation of current sheets either at the surface, or at 
the crust-core interface. The purely vertical drift acting within a 
spherical shell (the crust) will result inevitably in fast dissipation 
locally, wherever current sheets are formed. Moreover, when we 
consider NS models with a stratified electron number density, 
the gradient in n e through the NS crust may additionally support 
the creation of current sheets (mainly close to the crust-core in- 
terface). 

But there is yet another important problem related to the 
change from parabolic to hyperbolic character of the induction 
equation for purely toroidal fields. It is well known that the so- 
lutions of the inviscid Burgers equation develop discontinuities 
or, if a small viscous term is present, steep gradients. This is 
probably the reason why spectral methods have always failed 
when trying to solve numerically the induction equation for 
large magnetization parameters. The system evolves naturally to 
form a discontinuity or sharp gradient in B^ on a characteristic 
timescale which can be shown to be « th 3 ii (the Hall timescale 
is defined below in Eq. d30b ) 



dB u 



^"shock 



,2nen e R dz 



Rfr IdB, 



2f \ dz 



(26) 



which is of the order of * 10 6 - 10 7 yr for B = 10 13 G and 
typical lengthscale of 1 km. However, it can be much shorter 
for stronger fields and small scale structures generated by the 
Hall cascade. A similar timescale characterizes the typical travel 
time of a magnetic perturbation to reach either the surface or the 
core-crust boundary, 



5, 



^travel 



\2.nen e R d 



Table 1. Description of initial models, differing in the initial field 
strength and relative sign of the poloidal and toroidal compo- 
nents. All initial poloidal fields are purely dipolar (n = 1). The 
magnetic field is in units of 10 13 G. 



Model 




"tor 


Multipole 




(10 13 G) 


(10 13 G) 


(toroidal) 


A 


10 


-100 


n = 2 


B 


10 


100 


n = 2 


C 


10 


-100 


n = 1 


D 


1 


-10 


n = 2 


E 


20 


-200 


n = 2 


F 


50 


-50 


n = 1 


G 


100 


-100 


n = 2 



(27) 



where d denotes the distance between the initial position of the 
perturbation in the crust and the surface or the crust-core in- 
terface. The formation of a shock by the compression of the 
field against the inner or outer boundaries will cause unavoid- 
able numerical instabilities and/or the Gibb's phenomenon if one 
tries to solve the problem by means of spectral methods. Note 
that, locally and temporarily, the three characteristic timescales 
(i"Haib T s hock, and Ttravei) may coincide or differ by orders of mag- 
nitude, because they are proportional to the second, first, and 
zeroth derivative of the field strength, respectively. These char- 
acteristic timescales become different once small-scale structure 
or current sheets appear, creating small length scales in addition 
to the larger crust size. This makes any reasoning based on those 
timescales questionable; only a serious numerical study of the 
Hall drift can yield an idea about its effects on the crustal mag- 
netic field evolution. 

For the simple case of a Burgers-like equation with constant 
n e , (see Eq. (|25T>), we have checked that upwind methods, specif- 
ically designed to deal with hyperbolic equations, work very 
well in regimes «bt — > oo, in which a spectral method fails. 
Thus, the numerical problems observed with spectral codes for 
o>bt ^ 100 are most likely caused by an intrinsic limitation of 
the numerical approach. In the general case, having both poloidal 
and toroidal field components in a stratified spherical shell, the 
field evolution is not so simple. Then, the equations are strongly 
coupled, they have both a parabolic and a hyperbolic part and it 
is difficult to guess what is the best strategy to solve them. It is 
not the scope of this paper to give the final answer about the best 
numerical technique. We intend to point out in this section that 
there are deep unavoidable reasons that lead to unsurmountable 
problems in many cases. For the rest of the paper, we will fo- 
cus on simulations with realistic NS models that can be handled 
by our hybrid method. We take advantage of a seldom occasion: 
the more realistic model causes less numerical problems than 
the constant density model. The reason is that current sheets are 
faster smoothed out in the shell layers just below the surface, be- 
cause their electric conductivity is orders of magnitudes smaller 
than in the inner crust. In addition, in the long run, the toroidal 
field seems to find a quasi-equilibrium configuration in which 
the effect of the gradient of conductivity is counterbalancing the 
Hall drift. Such an equilibrium can not be established for non- 
stratified crustal models. 



4. The NS model and the initial conditions. 

Our aim is to study the global evolution of the magnetic field in 
isolated NSs. The Hall drift occurs both in the fluid core and in 
the solidified crust. While the effect of the Hall drift in the core 
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10 kyr 50 kyr 500 kyr 




Fig. 2. Top: Poloidal magnetic field lines for three different evolutionary stages of the magnetic field configuration initially described 
by model A. (f = 10, 50, and 500 kyr). Bottom: contours of constant toroidal field strength (B^). The crustal shell has been stretched 
a factor of 4 for clarity. 



is less obvious and may proceed on a timescale of the order of 
the Hubble time (insert L5 = 10 and p « 10p nuc into Eq. (61) of 
Goldreich & Reisenegger (1992)), its effect can be crucial in the 
crustal field. 

In order to build the background NS model we have used a 
modern Skyrme-type equation of state (EOS) at zero tempera- 
ture describing both, the NS crust a nd the liquid core, based o n 
the effective nuclear interaction SLy (Douc hrn & Haensel 1991b . 
The low density EOS (below the neutron drip point) employed 
is that of Baym et al. (1971). With this EOS, we have built 
a NS model having a radius of about 11.7 km and a mass of 
1.28 M Q . The central density is 8.83xl0 14 g/cm 3 and its crust 
(from p a 10 10 g/cm 3 to ss 10 14 g/cm 3 ) extends from 10.7 
to 11.6 km. In Fig. Q] we show radial profiles of the electric 
conductivity <x (solid lines) and <jl>btIB\t, (dashed lines), where 
Z?i3 is B in units of 10 13 G, for three different temperatures. 
We have assumed a constant impurity concentration parameter 
Q = Himp(Z im p - Z) 2 /n ; of Q 



10 2 . Here, n lmp and n-, are the 



impurity and ion particle densities, respectively, and Z lmp is the 
charge number of the impurities. We have compared simulations 
with Q = 10 -2 and Q - \0 A without finding significant differ- 
ence s. A high im purity content could lead to even faster dissipa- 
tion (I Jones 20041) . The figure illustrates the fact that the electric 
conductivity varies by 3 - 4 orders of magnitude within the crust 
and depends strongly on the temperature. The magnetization pa- 
rameter scales linearly with B. For the fiducial field of 10 13 G, 
it is of the order of unity for a temperature of 10 9 K but can be- 
come as large as 1000 as the star cools down. For magnetar field 
strength which both at the surface (radial poloidal component) 



and within the crust (meridional poloidal and toroidal compo- 
nents) may well exceed B13 = 10, the magnetization parameter 
can reach locally values in excess of 10000. 

From Fig. [TJone can also read the relevant timescales of the 
problem. The Ohmic timescale is 



47TCr^ 2 



T"Ohn 



(28) 



where A is the typical magnetic field length-scale. Inserting some 
typical numbers, we get 



TOhm =4.4 



,10 2 V 

On the other hand, the Hall timescale is 



1"Hall 



cB 



(29) 



(30) 



The ratio of the Ohmic to Hall timescale is simply given by u>bT. 
From Fig.Q]we can infer that, if the temperature of a NS's crust 
varies between 10 9 and 10 8 K during the first million years of 
its life, the average Ohmic timescale in the crust is « 1 Myr. 
For magnetized NSs, the ratio of the Ohmic to Hall timescale is 
approximately 



TOhm 
THall 



(1 - 10) Xfi! 



(31) 



Note, however, that such averaged timescales are of very re- 
stricted use to characterize the field evolution in NS crusts, since 
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Fig. 3. Same as Figure[2]but for model B. 



both the density and the magnetic field vary over many orders 
of magnitude there. There is no doubt, then, that in any circum- 
stances there is magnetic field decay in the crust of young NSs. 
A different issue is whether or not this effect is observable when 
studying populations of older NSs. Because of the strong tem- 
perature dependence of the electrical conductivity, when a NS's 
crust cools below 10 7 K, typically one million years after birth, 
the Ohmic dissipation time increases significantly, and no rapid 
field decay can be expected after that age. 

For this reason, in this paper we want to focus on the initial 
evolution of magnetic fields in relatively young NS. In order to 
mimic the cooling process of a NS, we have started with a crust 
at T = 10 9 K, a typical value after formation of the crust, at 
most within hours after birth, and we force the temperature of 
the isothermal crust to vary according to 

T(t) = 10 9 (1 + 10 6 r 6 r I/6 (32) 

where is the NS age in 10 6 years. This approximation is valid 
during the neutrino cooling era if only modified URCA pro- 
cesses are operating (Page et al. 2006). It is not the purpose of 
this paper to discuss how the thermal and magnetic field evolu- 
tion are coupled in detail, but this simple approximation is suffi- 
cient to capture the main effect: as the NS's crust cools (from 10 9 
to 10 8 K in about 1 Myr), via the decrease of the electron relax- 
ation time (r), the Hall term becomes more and more important. 
When a significant part of the crustal field is dissipated and/or it 
has approached a force free configuration, its decay continues on 
a much longer Ohmic timescale. We have not included effects of 
temperature anisotropies within the crust, although they may be 
important. This will be addressed in detail in future works. 



The other microphysical input needed for performing the 
simulations is the electrical conductivity and the magnetization 
parameter. Since, as the temperature drops and the magnetic field 
evolves, the magnetization parameter and the electrical conduc- 
tivities vary, at each time step we recompute the value of the re- 
laxation time and the electrical conductivities by using the public 
code developed by A. Potekhin (1999)0. Despite it causes some 
numerical efforts, we have chosen to be as realistic as possible 
and use this state-of-the-art microphysical ingredients to account 
for the effect of composition stratifications. 

4.1. Initial magnetic field configuration. 

Little is known with certainty about the initial magnetic field 
strength and structure. Recently, Braithwaite (2006) and Geppert 
& Rheinhardt (2006) have shown that sufficiently rapid rota- 
tion can stabilize dipolar toroidal and poloidal fields of magnetar 
strength against MHD instabilities occurring immediately after 
completion of the proto-NS phase. Only preliminary results are 
available about the stability of magnetic field structures consist- 
ing of a poloidal-toroidal mixture. Our conclusion about the ori- 
gin of strong toroidal crustal fields that probably exist in mag- 
netars and quite certainly exist in thermally emitting NSs with 
a highly non-uniform surface temperature distribution (Geppert 
et al. 2004, 2006; Perez-Azorfn et al. 2006a, 2006b) is, there- 
fore, that they have probably been generated at birth. Later, they 
were stabilized against Tayler instabilities, and have been frozen 
into the solid crust when it forms. Alternatively, the presence of 
sufficiently strong temperature gradients both close to the crust- 

2 www. iof f e . rssi . ru/astro/conduct/condmag . html 
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Fig. 4. Same as Figure|2]but for model C. 



core interface and in the degenerate surface layers (Wiebicke & 
Geppert 1995) are able to convert thermal into magnetic energy 
very effectively via a thermoelectric instability. In this case, the 
magnetic energy will be predominantly stored in toroidal field 
structures. 

We have considered quite different initial structures: purely 
poloidal and purely toroidal ones, with initial dipolar and 
quadrupolar modes, and several mixed initial fields where the 
ratio between the energies stored in the toroidal and poloidal 
field parts is varied. For the poloidal field, our initial configu- 
rations are chosen to be free Ohmic decaying modes (strictly 
speaking, pure decay modes for constant conductivity profiles) 
satisfying the boundary conditions. For a given angular multi- 
pole (n), the radial eigenfunctions can be written as linear com- 
bination of spherical Bessel functions of the first (j„) and second 
(n„) kind 



<E>„0) = a r j„(jur) + b r n„(/zr) 



(33) 



where \x is the radial wavenumber. The boundary conditions fix 
the wavenumber fi and the ratio b/a. For the models presented 
in this paper, the overall normalization factor has been chosen 
to fix the radial component at the pole of the dipolar (n — 1) 
component to the values of B po i listed in Table [1] 
The initial toroidal field is given by 



TO = c[(r-/?,)(r-fl)] 2 



(34) 



where the constant c is determined by the maximum value of 
B v . In Table 1 we summarize the models employed in this work, 
differing in the initial values and relative sign of the poloidal and 
toroidal components. 



5. Results and discussion. 

We present now the results of our numerical simulations of the 
field evolution in the crust and its dependence on the initial mag- 
netic field structure and strength. We will restrict this presenta- 
tion to the selection of initial models described in Table [1] al- 
though we have performed a number of different simulations 
with a variety of initial configurations. In particular, we have re- 
produced some of the toy models (cons tant & and t) found in the 
literature (Hol lerbach & Riidiger 2002 ) to test our code. We saw 
a good qualitative agreement with only minor quantitative dif- 
ferences. We have also performed simulations of the evolution 
of purely toroidal fields, which remain toroidal forever, for test- 
ing purposes. For conciseness, we will discuss now only realistic 
models of magnetized NSs. 

In Fig. |2] [3] and |4] we present the evolution of models A, B 
and C. For three representative ages, we show both the poloidal 
magnetic field lines (top panels) and the contours of constant 
toroidal field strength (bottom panels). The crustal shell has been 
stretched by a factor of 4 for clarity. Models A and B differ only 
in the relative sign of its toroidal component which is for both 
models initially quadrupolar (n = 2), while model C has an ini- 
tially dipolar (n = 1) toroidal field with the same radial profile 
as model A. In all cases the models have the same maximum 
toroidal field strength. The different sign of the toroidal field 
affects the direction of the Hall drift (toward north or south). 
The other important effect is due to the gradient of the electric 
conductivity, which always causes the tendency to displace field 
lines towards the inner crust. 

For Model A, we see in the central panel how the Hall drift 
compresses magnetic field lines near the equator, while the re- 
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Fig. 5. Radial profiles of <t>i and (in arbitrary units). These 
quantities are proportional to the poloidal (dipolar) and toroidal 
(quadrupolar) components, respectively (see Eq. |6j. The upper 
panel shows the initial model, while the central and bottom panel 
correspond to t = 20 kyr and t = 500 kyr, respectively. The dif- 
ferent line styles stand for models A (solid line), B (dotted line), 
and E (dashed line). Notice the different normalization of rT^ in 
the upper and middle panel (a factor of 10 and 5, respectively). 



suits for Model B show the opposite tendency. Model C, due 
to its dominant initial dipolar structure, shows a global dis- 
placement of the magnetic field towards the south pole, that 
results in a fast dissipation in the high resistivity surface re- 
gion near the south pole. During that phase, the interplay of the 
poloidal and toroidal field modes is very intense. The magnetic 
field is dragged and twisted, thereby creating current sheets. At 
these sites of very efficient dissipation the toroidal field weakens 
rapidly and, after becoming about equally strong as the poloidal 
field, the latter tends to return into its original position. This 
sequence of twisting and stretching and subsequent release of 
the twisted field characterizes this epoch by its oscillatory be- 
haviour. However, looking at the right panels, we see that after 
about half a million years, all poloidal fields have a very similar 
appearance. The strongly dissipated, and now weaker, toroidal 
components, become more regular than in the Hall phase, and 
they are dominated by the n = 2 or n = 3 modes. 

The main conclusion from these results is that the initial 
magnetic field configuration largely determines the early evolu- 
tion (< 10 s yr). On the long run, however, there seems to exist a 
clear tendency to establish a more stable configuration, consist- 
ing of a dipolar poloidal component in combination with higher 
order toroidal field modes. This fact can be better understood 
when looking at the evolution of the radial profiles of the dom- 
inant modes, as shown in Fig. At early times, say t — 20 
kyr, the interplay between the toroidal and poloidal components 
results in a different evolution of <J>i, although the shape is qual- 
itatively similar. After half a million years, the poloidal field is 
clearly dominant, and very similar for all models. The toroidal 
field is weaker (in real units, about a factor of 2-3 smaller than 
the poloidal field) and concentrated towards the inner crust. It 
is interesting to notice that also ¥2 is rapidly driven towards a 
similar shape in all cases, indicating that after the fast initial 
transient the gradient of the conductivity determines a sort of 
quasi- equilibrium field. 

It is also interesting to look at the evolution of the total mag- 
netic energy (Fig. |6]l to understand the different evolutionary 
stages that the NS passes through. Firstly, it is remarkable how 
all models have a very similar evolution, and how they converge 
at late times towards a similar exponential decay. For all mod- 
els, the magnetic energy has been dissipated by a factor 10 to 
20 after 0.5 Myr. Secondly, we can clearly differentiate between 
the initial fast decay, in which the field can dissipate one half of 
its energy in only ~ 10 4 yr. This is caused by the much smaller 
initial conductivity at higher temperatures combined with the ef- 
fects of the Hall drift. During the subsequent stage, the Ohmic 
dissipation timescale (about half a million years) is almost iden- 
tical for all models. Obviously, the first stage is of great impor- 
tance for models with initially stronger fields, i.e. it will play an 
important role for NSs born as magnetars. This seems to be the 
only period, before the photon cooling era, when the Hall drift 
will qualitatively affect the field evolution. In general, the mod- 
els in which the initial ratio of toroidal to poloidal field is small 
(models F,G) decay slower than those with large ratios (mod- 
els A,B,C,E). Among the latter, the stronger the initial field, the 
more efficient is the dissipation during the Hall phase. 



3 The black/white figures presented here are snapshots and yield 
only a vague impression of the evolution of the magnetic field with 
Hall drift and Ohmic diffusion, everything coupled with the cooling of 
the crust. We advice the reader to look at the movies, available under 
http://www.dfa.ua.es/UNS06/movies.html , where the evolution 
is visualized in much more detail. 
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Fig. 6. Evolution of the total energy stored in the magnetic field 
formalized to the initial value, for models A (solid), B (dotted), 
C (short dashes), E (dot-dash), F (triple dot-dash), and G(long 
dashes). The thin line shows, for comparison, an exponential de- 
cay with a timescale of « 0.4 Myr. 
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Fig. 7. Evolution of the total energy stored in the magnetic field 
for models A (solid) and G (dashed) compared with the same 
initial configurations evolving only by Ohmic decay (thin lines). 
The difference between thick and thin lines is due to the nonlin- 
ear Hall terms, that result in a faster initial dissipation. 



The Ohmic dissipation rate appears to be faster than what 
has been usually discussed in the literature. The reason is that 
we are simulating for the first time the evolution of NS models 
with a realistic profile of electron density, composition, conduc- 
tivity, etc., including the coupling between different field modes 
through the non-linear Hall term. The effect of the Hall drift 
is twofold: it causes the transfer of magnetic energy to smaller 
scale modes that decay faster and, in some cases, the displace- 
ment of the magnetic field to regions of higher resistivity where 
it is rapidly dissipated. 

In Fig. [7] we compare the evolution of the total magnetic en- 
ergy in models A and G (thick dashes and solid lines, respec- 
tively) with the hypothetical evolution of the same initial config- 
uration without the Hall drift (by setting f = 0). After the initial 
phase (about 50 kyr) during which the models that include the 
Hall drift dissipate energy about 10 - 20% faster, the rate of en- 
ergy dissipation becomes approximately the same. This makes 
evident that the initial departure of the exponential decay is par- 
tially caused by Hall drift, but also by the relatively rapid cooling 
of the crust that results in a significant time-dependence of the 
conductivity and of the magnetization parameter. Other cooling 
models, for example fast cooling due to the activation of direct 
URCA processes, may result in an initially different evolution. 
However, we think that on a secular timescale of about 0. 1 Myr 
the system readjusts itself and the dissipation of magnetic field 
is mostly controlled by the Ohmic decay. In a future work we 
plan to couple a multidimensional cooling code with the mag- 
netic field evolution to study different cooling scenarios. 

To finish the discussion in this section, we present in Fig. [8] 
the power spectra for models A, B, and C at three different mo- 
ments of the evolution. In the left panel, which corresponds to 
t — 1 kyr, we see how the Hall cascade is filling very quickly 
the shorter wavelength modes (initially only n = 1 and n = 2 
modes existed). At t = 10 kyr, which corresponds roughly to the 
Hall timescale, the Hall cascade and perhaps the Hall instability 
have filled out all large wavenumber modes and approximately 
saturates following a n~ 2 power law. This situation is kept for 
about 50 kyr (not shown in the figure), until the field has been 
dissipated by a significant amount, and the Hall term begun to 
lose its importance with respect to the regular Ohmic dissipa- 



tion term. After half a million years the power spectrum is much 
steeper (oc n~ 6 ), an indication that the Hall drift has lost its influ- 
ence on the crustal field evolution. 

5.1. Spin down evolution. 

Since the evolving crustal field may also have effects on the ro- 
tational evolution of the NS, we study the differences of the spin 
down behaviour between different models. In Fig. [9] we show 
the evolution of the period (top), period derivative (middle) and 
spin-down age (f s( j = P/(2P)) (bottom) in four selected models 
with different initial dipolar field strengths. The line-styles cor- 
respond to models A (solid lines), D (dash-dotted line), E (dots), 
and G (dashes). A typical NS born as a magnetar would corre- 
spond to model G or, perhaps, E, while a typical pulsar such as 
Geminga would show the behaviour of model D. 

Notice that a minimum initial poloidal field of about 2 x 10 14 
G is required to explain the large rotation periods (5 - 10 s) of 
thermally emitting isolated NSs as the Magnificent Seven. This 
may be another indication towards a common evolutionary path 
in which some INSs were born as magnetars and their magnetic 
field has decreased one order of magnitude during their life- 
times. The evidence for crustal field decay presented here has 
also implications for estimates of the ages of old NSs. Pulsar 
spin down is thought to generally follow the prediction of the 
vacuum dipole model, which gives v oc B 2 v 3 , where v is the spin 
rate. If the birth spin rate far exceeds the present spin rate and B 
is constant, the age in this model is given by f si j. This expression 
is used as the standard estimate of a pulsar's age. If the field is 
decaying according to our simulations, however, the relationship 
between the true age and f s d given by the dipole model is shown 
in the bottom panel of Fig. [9] It shows that f S( j seriously overes- 
timates the age for NS older than 10 5 years. This effect helps to 
reconcile the observed discrepancy between the spin-down ages 
and independent measures of the ages of some isolated NSs. 

Another interesting feature is the rapid variation of the dipo- 
lar poloidal component during the Hall epoch. Of course, the 
total magnetic energy always decreases, but the n — 1 poloidal 
component may exhibit an oscillatory behavior. This is a result 
of the efficient energy transfer between the large scale toroidal 




and the dipolar poloidal field modes, i.e. an genuine effect of the 
Hall drift. Assuming that only the n — 1 component contributes 
to spin down the star, the oscillations of the dipolar poloidal sur- 
face field affect the braking index n = vv/v 2 , as shown in Fig.fTUl 
During the first ~ 10 5 years of their lives, sufficiently magnetized 
NSs will show a quite erratic variation of the braking index, that 
may reach any value from 1 .5 to 4. The oscillations of the brak- 
ing index continue even until ~ 5 x 10 5 years, when the initial 
toroidal field is stronger, as seen for model E which has an initial 
maximum toroidal field of 2 x 10 15 G. At late times, once Ohmic 
dissipation is controlling the evolution and the dipolar magnetic 
field decreases steadily, the braking index gradually increases 
and might reach very high values. This process will be inverted 
after the crustal field is dissipated almost completely and the ro- 
tational evolution of the NS is determined by the much longer 
timescale of the core field expulsion. In that period, the brak- 
ing index will approach its canonical (dipolar) value of 3 again. 
Our results for the braking index evolution coincide very well 
with the observed indices for middle aged pulsars (see Geppert 
& Rheinhardt g002). It can be shown that the braking index sat- 
isfies the relation 

n = 3 -2— k 3 +4— , (35) 
BP t b 

where we have introduced a typical magnetic field decay time 
t b = B/Bto derive the latter equality. The thin solid line in Fig. 
[10]corresponds to this approximation with t# = 0.8 Myr. Notice 
that Fig. [7] shows the evolution of the total magnetic energy (oc 
B 2 ), so that the decay time for the magnetic field is twice the 
value quoted in Fig. [7] 

It should be noted that after « 1 Myr neutrino emission no 
longer controls the thermal evolution of a NSs, and our results do 
not apply. Moreover, the low temperatures reached at that time 
would increase significantly the conductivity (see Fig. [TJ. Then, 
the purely Ohmic decay will proceed much slower, but also the 
magnetization parameter will increase, which may result in an- 
other Hall stage during the photon cooling era. In future work, 
we plan to extend our study to longer times with consistent tem- 
perature evolution. In addition, one always must keep in mind 
that the magnetic field component supported by currents and in 
the superconducting NS core would be dominant after the crustal 



field has been dissipated. This, however, is a completely differ- 
ent scenario that cannot be analyzed with our presently available 
tools. 

5.2. The Hall instability ? 

We have seen in our numerical studies some hints for the oc- 
currence of the Hall instability, i.e., for the rapid, non-local in 
the momentum space, energy transfer from the initial large scale 
field modes into much smaller ones. Related spectral features 
can be seen e.g. in Fig. [8] where local maxima appear. However, 
at the present level of energy conservation (99.9%) and our lim- 
iting angular resolution (n max = 120, although most of our runs 
were performed with « max = 50) we can not clearly distin- 
guish the truncation effects from the onset of the Hall instability. 
Therefore, we decided to postpone the study of the Hall instabil- 
ity in the NS crust to future work. 



6. Conclusions 

The above presented results show that magnetic fields main- 
tained by currents circulating in the crust of NSs are strongly re- 
arranged and do decay significantly during the first million years 
of a NS's life. In addition to purely Ohmic decay, which is faster 
in the first thousands of years when the electrical conductivity 
is relatively large due to the high crustal temperature, we find 
that the Hall drift may contribute noticeably to accelerate the 
dissipation of magnetic fields. For typical field strengths of 10 14 
G, we observe a Hall drift dominated stage followed by purely 
Ohmic decay proceeding in a timescale of the order of 1 Myr. 
Depending on the strength and structure of the initial magnetic 
field, this Hall phase lasts a few 10 3 — 10 4 years and it is char- 
acterized by an intense exchange of magnetic energy between 
the poloidal and toroidal components of the field and by the re- 
distribution of magnetic field energy between different scales. 
It can be expected that such rearrangements and the relatively 
rapid field decay have observational consequences, as can be ob- 
served in magnetars. Of course, whether this first phase plays an 
important role and how distinctive it is depends on the initial 
field strength and structure and how fast the NS cools. If a NS 
begins its life as a magnetar its external dipolar field is > 10 14 



12 



J.A. Pons & U. Geppert: Magnetic field dissipation in neutron star crusts 



10.0 r 



1.0 r 



0.1 r 




0.01 0.10 

t (Myr) 



1.00 



10 -io 




0.01 0.10 
t (Myr) 



1.00 



1.000 



0.100 r 



0.010 r 



0.001 





,y ' "y- 

// / '■ 


■ J? 




yf ' 




fi 


i 



0.01 



0.10 



1.00 



t (Myr) 



Fig. 9. Evolution of the period (top), period derivative (middle) 
and spin-down age (bottom) in four selected models with dif- 
ferent initial dipolar field strengths. The line-styles correspond 
to models A (solid lines), D (dash-dotted line), E (dots), and G 
(dashes). A typical NS born as a magnetar would correspond to 
model G, while a typical pulsar such as Geminga would show 
the behaviour of model D. 




Fig. 10. Evolution of the braking index n - vv/v 2 for the same 
models as in Fig. [9] 



G. Within the crust, however, the internal magnetic field may lo- 
cally exceed that value by about one order of magnitude. Thus, 
even the expected initial high crustal temperatures of ^ 10 9 K, 
which cause a relative small electron relaxation time (r), cannot 
avoid that wgt » 1 in a large fraction of the crust volume. 

The toroidal part of the field is specially affected by the Hall 
drift. There are two main effects acting upon the toroidal fields: 
it is globally displaced toward the inner crust because of the al- 
most everywhere negative conductivity gradient and, depending 
on the relative sign with respect to the poloidal component, it 
tends to move vertically toward one or the other magnetic pole. 
After the Hall stage, during which the toroidal field is strongly 
rearranged and dissipated, the long term evolution seems to se- 
lect, generally, a predominantly quadrupolar/octupolar structure 
concentrated in the inner crust and with tendency to be stronger 
close to the poles. This multipolar structure will determine the 
surface temperature distribution of middle aged NSs, that could 
be more complex than previously thought. Such complex field 
structures and the local deposition of energy by Joule heating fa- 
vor surface temperature distributions characterized not only by 
two hot polar caps, but for example by a hot equatorial belt, 
as has been probably seen in RX J0720.4-3125 (Haberl et al. 
2006, Perez-Azorin et al. 2006b) or can be inferred from the 
light curve of RBS 1223 (Schwope et al. 2005). 

Though the toroidal part of the crustal field undergoes a 
spectacular dissipation and rearrangement, the coupling between 
both parts also affects the evolution of the poloidal field, which 
is responsible for the spin down of NSs. Thus, the complex in- 
teraction of Hall drift and Ohmic dissipation is also reflected in 
the temporal behaviour of the braking index and, in principle, 
accessible for observations. Our models show that the braking 
index of young (;S 10 5 years), magnetized NSs exhibits a quite 
wild variation, and can reach any value from 1 .5 to 4. After the 
Hall stage, when Ohmic dissipation controls the evolution and 
the dipolar magnetic field decreases steadily, the braking index 
gradually increases and might reach very high values. This pro- 
cess will cease after the crustal field is dissipated almost com- 
pletely, or when the conductivity increases as the NS cools down. 
During the photon cooling era, the braking index will approach 
its canonical (dipolar) value of 3 again. 

If the initial magnetic field is too weak (^ 10 12 G) for the 
Hall stage to be relevant, the evolution will proceed according to 
purely Ohmic field decay. The very existence of magnetars and 
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of their presumable descendants, thermally emitting NSs with Appendix A: Induction equation in cylindrical 

B « 10 13 G, suggests that the fraction of NSs born with large coordinates, 
magnetic fields may be larger than expected. 

In this appendix we use a different set of variables to derive the 
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de Ciencia y Tecnologi'a grant AYA 2004-08067-C03-02. JAP is supported by a r cos ff) an( J use t he notation of Hollerbach & Rudiger (2002). 

Ramon y Cajal contract from the Spanish MEC. Alternatively to our decomposition in Eq. ©, one can in axial 

symmetry simply work with the ^-components of the vector po- 
tential and of the magnetic field, 

B = B pol + B t01 = V x (A^) + B v e^ . (A. 1) 

In cylindrical coordinates, we can write explicitly the different 
components of the magnetic field and the current density: 

dA„ 1 d(RAJ 
B = ~~dF eR + R~d^ £: + B ^ ' (A ' 2) 

^ „ dB « id(RB^) 
— J = VxB = ~^re R + ~^^e z + , (A.3) 
c oz R oR 

where we have introduced the notation 



J\ = -V 2 A^ + -± . (A.4) 
The induction equation in terms of this variables reads: 

= [( V X *tar) X flpol] - 4- 

Ot cr cr 

dB v It \ 

= «(p-Vxl-[(Vx Bpoi) x B pol + (V x B tor ) x B tor ] I 



+ Vx 



IvxB^ 

cr 



(A.5) 



The term appearing in the r.h.s. of the first equation can be writ- 
ten as follows: 

e v ■ [(V x B tor ) x B^] = i [V(RAp) x V(RBj] ■ (A.6) 
The two terms inside the curl of the second equation are 

(Vx B tor ) X B te = —£V(RB V ) 
R 

(VxB tor )xB pol = ^V(RA^) (A.7) 

K 

Then, taking the curl and after some algebra, the two terms ap- 
pearing in the equation for B^ can be written as follows: 

■ V x [(V x B tor ) x Bte] = [V(/i8 f ) x V(B f /fi)]^ 

2B^ dBy 
~ ~~R~~dz 

c,-Vx[(Vx B pol ) x Bpoi] = - [ViRAp) x V( , (A.8) 
where the first one is used in the derivation of Eq. d25l l. 
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